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Abstract 

For networks of coupled dynamical systems we characterize admissible functions, 
that is, functions whose gradient is an admissible vector field. The schematic represen¬ 
tation of a gradient network dynamical system is of an undirected cell graph, and we 
use tools from graph theory to deduce the general form of such functions, relating it to 
the topological structure of the graph defining the network. The coupling of pairs of 
dynamical systems cells is represented by edges of the graph, and from spectral graph 
theory we detect the existence and nature of equilibria of the gradient system from the 
critical points of the coupling function. In particular, we study fully synchronous and 2- 
state patterns of equilibria on regular graphs.These are two special types of equilibrium 
configurations for gradient networks. We also investigate equilibrium configurations of 
S 1 -invariant admissible functions on a ring of cells. 


1 Introduction 

Networks of coupled dynamical systems frequently demonstrate phenomena, such as syn¬ 
chrony, phase relations, resonances and non-generic bifurcations, that are typically associ¬ 
ated with group equivariance, even when the networks are not invariant under any non-trivial 
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group actions. They appear to have some form of ‘hidden symmetry’. These phenomena 
have been explored from a number of different points of view, e.g. using groupoid formalism, 
as in [ 9 ] and m or, more generally, category theory, as in [ 6 ], and semigroup formalism, as 
in D2I- The algebraic formalism of symmetry groupoids of networks has led to a successful 
way to search for patterns of synchrony on networks. In general, elements of the groupoid 
can be thought of as a set of local symmetries of the network, relating sets of cells in the 
network in a way that the system is invariant under the action of these symmetries. This 
has been first established in |19j and explored since in a great number of works; for example, 
[21171122], among many others. 

In this paper we focus on one particular class of systems, namely gradient systems. 
Synchrony manifests itself in a network as configurations of the coupled systems that behave 
setwise identically. For gradient systems there are also configurations given by phase relations 
between cells, which appear naturally in a class of vector fields generated by functions with 
an extra S 1 -invariance. 

We consider a network of a finite set of dynamical systems cells coupled together in a 
manner given schematically by a connected undirected graph Q = (V,S) without multiple 
edges, whose vertices represent the cells and edges correspond to couplings. We shall call 
Q a cell graph. The set V = {ui,..., v n } denotes the set of vertices and £ the set of edges 
of Q, which are identified as 2-element subsets of V although we shall denote an element of 
£ as a pair (ui,u 2 ). We assume that a cell graph may contain all internal edges (or loops), 
in which case VCf. We shall also denote by I(v) the input set of v G V, that is, the set 
of vertices u G V, u ^ v , such that (u,v) G £, and d(v) shall denote the degree of v, the 
number of vertices in I(v). Notice that for our purposes we shall not consider loops to be 
in the input sets. A smooth manifold P v is assigned to each vertex v G V, so that the total 
configuration space is P = JX e y Pv- We denote this network by Af = (Q, P ). 

A vector held g — (g ±,..., g n ) on P is called an admissible vector held if it is consistent 
with the network structure, that is, if it is equivariant under the action of the network 
symmetry groupoid (Definition 4.1 in [TUj 1. In particular, if x v denotes coordinates on P v , 
each component g v must depend only on the variables x u for which (u, v) G £. More precisely, 
for I(v) = {ui,... ,u d( „)}, 

9v(%) = £iv{x v , x vi ,..., ) • (1 • 1) 

In addition, equality constraints are imposed between components corresponding to cells u 
and v whose input sets I(u) and I(v) are isomorphic by an element of the groupoid, that is, 

9u 9v (1-2) 

General network dynamical systems are usually defined for directed graphs. For a network 
of smooth dynamical systems to be given by the (negative) gradient of a smooth function 
/ : P —> R, 

x = —V/(x), (1.3) 

it is necessary that if (u,v) is an edge then so is (v,u), and hence the directed graph is 
equivalent to an undirected graph. We assume that all edges are identical and all vertices 
are of the same type. In particular, all components P v of P are the same. Since V/ is an 
admissible vector held, then its components satisfy (11.11) and (11.21) . Since we are assuming 
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that all cells are identical, then (II.2ft holds precisely if d[u) = d(v), and, for each v G V, 


df_ 

dx v 


(x) 


fd(v) 5 % 


u 1 7 


j Xu d(v))l 


(1.4) 


where the over-bar indicates invariance by permutation of these variables. 

Simple examples show that not all admissible vector fields are gradient. In this paper 
we restrict to the class of admissible vector fields that are of gradient type. Onr main 
result, Theorem 12.41 characterizes the admissible functions, namely, smooth functions whose 
gradients are admissible vector fields. As we shall see, for any network these functions 
are decomposed as a sum of components that depend on each cell individually - the self¬ 
connection functions - and those that depend on the way the cells are coupled together - the 
coupling functions. 

This work also addresses the analysis of critical points of admissible functions / or, 
equivalently, the analysis of equilibria of (1 1.3 1) . As suggested by the general form of these 
functions, there is a direct relationship between critical points of the associated coupling 
function and critical points of the admissible function that are either totally synchronous 
(all cells assume the same value) or given by 2-colour patterns (when each cell assumes one 
out of two possible distinct values). In fact, this makes these two configurations on networks 
special in the class of gradient systems. In this paper we investigate these two particular 
types of critical points, regarding existence and nature, and shall understand how these are 
related to both the critical points of the coupling functions and the architecture of the graph 
of couplings. This study does not generalize for patterns with more than two colours. For 
these cases a different general approach can be applied, and this is done in [I]. In the presence 
of extra symmetries, however, other types of critical points of admissible functions become 
expected for some classes of graphs. These can appear with a variety of configurations, with 
no nontrivial synchrony among cells for example, but still with a phase relation. Specifically, 
this is the case when the admissible function is invariant under the circle group S 1 . 

There are similarities between the network gradient dynamical systems that we consider 
and Ising-Potts model, Kuramoto model, antiferromagnetic XY model (also called AFXY 
model), neural networks and other interaction systems. For example, applications of our 
results given in Section [4] are related to the results of [8], where the author applies the method 
of averaging to reduce models of discrete and a continuum of neural arrays to systems of 
phase equations, that is, to equations in which each cell is represented by a single variable 
lying on S 1 and interactions between two connected neural cells are periodic functions that 
depend only on the difference between their two phases. In such model, an edge (u , v ) G £ 
represents a synapse between the two neurons u and v. There is also a direct connection 
between our study of critical points of admissible functions given in Subsection 14.21 and 
several results about ground states in Kuramoto and AFXY models (see BQ E21 m IB]). 

This paper is organized as follows: In Section 0 we present the general form of admissible 
functions on networks, which is our main result. In Section [3] we study critical points 
of admissible functions, relating their existence and nature with the corresponding critical 
points of the coupling functions. A detailed analysis of fully synchronous critical points 
and 2-colour patterns of critical points is carried out for the special case of regular graphs. 
Section [4] presents a study of all possible critical configurations of a class of admissible 
functions under an extra SMnvariance. 
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2 Admissible functions 


For a given cell graph Q = (V, £), we characterize the functions on P = W v& \> P v whose 
gradients are admissible vector held for Af = ( Q,P ). We shall consider P v = k > 1, for 
all v G V, taking local coordinates if necessary. The variable on each cell v G V shall be 
denoted by x v = (x \,..., xf). 

Definition 2.1. For a cell network Af = (Q,P), a smooth function f : P —>■ M is an 
admissible function if its gradient V/ is an admissible vector field for Af. 


The next result is one of the main constraints imposed on admissible vector fields if they 
are of gradient type. 


Lemma 2.2. For a cell network Af = (Q, P ) with |V| >3, if a smooth function f : P —y M 
is an admissible function, then, unless coupling is all-to-all, we have, for any 1 < i,j, l < k, 


d 3 f 

dx^dxlzdx 1 ^ 


( 2 . 1 ) 


if V \, V‘ 2 , V's are distinct vertices in Q. 


Proof: We prove by induction on the number of vertices in Q. For simplicity the proof is 
carried out for k = 1, but the result holds equally for any higher dimension k > 1. For 
|V| = 3, if (12.ip does not hold, then for each i, i = 1,2,3, df/dx Vi is a nontrivial function 
of the other two variables, and so Q necessarily has its three vertices coupled one another. 
For n > 4, assume that if a network has n — 1 cells and three distinct cells for which (12.ip 
does not hold, then these are all coupled one another. Now, let |V| = n and suppose that 
there exist distinct Vi,V 2 ,v 3 G V such that (12.ip fails. Let « 6 V, u / V\ ,u 2 , t’ 3 . Since Q is 
connected, there exists w G V such that (u , w ) G £. By the induction assumption all vertices 
in Q, except possibly u, are all-to-all coupled. In particular w G I(vf), i = 1,2,3. Now, for 
i = 1,2,3, we use the form (II. 4p of df/dx Vi together with the hypothesis of nonvanishing 
third-order derivatives with respect to x Vl ,x V2 ,x V3 to conclude that 


dx vi dx V2 dx w 


# 0 . 


( 2 . 2 ) 


But by the induction hypothesis we also have 


df_ 

dx w 


In-l(x 


w 5 ; 




,X, 


;)• 


So for all i — 1,..., n — 2, 


d 2 

dx u dx w 



) # 0, 


and hence u G D y/(u), that is, vertices in Q are all-to-all coupled. If u = ig, for some 
i G {1, 2, 3}, we just rearrange the indices in (12 .2 p to get the same conclusion. g 

We set up the following notation: for each edge e = (tq, v-2 ) G £ representing the coupling 
of vertices vi,V 2 G V, we denote p(e),r(e) G {ui^}, where p(e) ^ r(e) if v\ ^ u 2 and write 
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the ordered pair (p(e),r(e)) to represent a directed edge, with ‘head’ p(e) and ‘tail’ r(e), 
corresponding to the (undirected) edge e. 

As a consequence of Lemma 12.21 an admissible function is of the form 

f( x ) = ^M x p(e), x r(e)) + ^ a v (v), (2.3) 

ee£ ugV 

for smooth functions /3 e : M. 2k —> M and a v : —> M. We assume without loss of generality 

that the functions flfs do not contain terms of just one variable, and that / vanishes at the 
origin. For e G £ and v G V, the functions /Vs in the first sum shall be called coupling 
functions, and the a„’s shall be called self-connection functions. 

Theorem 12.41 below characterizes admissible functions defined on networks of coupled 
cells. As we shall see, these are distinguished when Q is a bipartite graph. 

Remark 2.3. Recall that a graph Q is bipartite if its set of vertices V can be devided into 
two disjoint subsets Vi and V 2 such that every edge of Q connects a vertex in a subset to 
a vertex in the other. In this case, a 2-colouring can be defined on the graph Q3j) such 
that any two vertices in either Vi or V 2 receive the same colour. For graphs of many cells, 
algorithms can be useful to check whether a graph is or is not bipartite (see HU and [IB]). 

For any subgraph S of Q let us denote £5 its set of edges. Also, let S n denote the 
permutation group acting on (R fc ) n by permutation of variables and let Z 2 denote the order- 
2 permutation group acting on any subspace (R fc ) 2 of two variables in (R^)" - . 

Theorem 2.4. If Af = (G,P) is an all-to-all coupling network of n cells, then admissible 
functions are the S n -invariant functions. Otherwise, a function f : P —>■ K is an admissible 
function associated to J\f if, and only if, there exist smooth functions ad(v) : R fc —>■ R and 
/3 : R 2fc —> R such that one of the following holds: 

1. If Q is bipartite, then for the disjoint union V = Vi U V 2 we have 

f( x ) = 5Z Pf x p(e)i x r(e)) 4“ ^ ^ ®d(v)( x v) T ^ ^ r )d(y)i x v)i ('^■^) 

ee£, p(e)eVi uelA 

where i = 1 or 2. In addition, if there exist V\ G Vi andv 2 G V 2 such that d{vi) = d{v 2 ), 
then j3 is Z 2 -invariant and 

f( x ) = J2 /3(-£p(e) > x r(e) ) T ^ ^ d(v) j, x v) • (2.5) 

e&£ v£\I 

2. If Q is not bipartite, then f is of the form H2.5\) with /3 Z 2 -invariant. 

Proof: If the cells in Q are all coupled one another, then V/ is S,,-equivariant, so / is 
S n -invariant. 

It is easy to verify that functions of the forms (12.4ft and (12.5p are admissible functions. 
For the converse, since we are interested in the non all-to-all coupling case, we have n > 3. 
We present the proof for k = 1, since for k > 1 the proof adapts straightforwardly. In fact, 
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the arguments rely essentially on results from graph theory and the repeated use of (11.4)1 : if 
k > 1, then for x v = {x \,..., x*), (|1.4|) is 


df_ 

dx l v 


(x) 


(v) ) X U\ 1 ■ ■ ■ 1 ) J i 1 ) 


and we work out the partial derivatives on each component of x v to obtain the same conclu¬ 
sions as for k — 1. 

Assume that Q is bipartite. Considering the partition V = Vi U V 2 , we have a decompo¬ 
sition 


S = 


where K? denotes the star graph K l d ( v \ with centre vertex v G Vi and d(v) edges (v,u), for 
all u G I(v). We notice that the disjoint union above refers to non-repeating edges, and it 
can obviously be taken to run through elements in either subset of the partition, indistinctly. 
Using the decomposition above, we rewrite (12.31) as 


/ 0 ) = P*( X fte)iX T {e)) + 

ueVi e&S/cv 


Q( ? , ( X v ) 


E E & 

ueVi uei(v) 


veV 


From (11,4)1 . for each v G Vi with d(v) > 2, the partial derivative df/dx v is invariant by 
pairwise permutation of u i: Uj G I(v). Hence, 


d/3 UiV 

dx v 


\x„ 


x ui ) 


d/3 Ujl 

OX,, 


( x v , x ui ). 


Since coupling functions do not depend on each variable independently, the equality above 
yields /3 UiV = f3 UjV . Thus, 


f(x) 


^ ^ ^ ^ ^ ^ v ipEy) • 

ugVi 


If Q is the star graph, then Vi = {u}. In this case, by (jl .4(1 . 

d P v / \ 1 / / \ dp v / \ 1 / / \ 

, x Ui ) + Q! u . (x Ui ) = — {x v , x Uj ) + a u . (x Uj ), 

for any pair u t , Uj G I{v). Set x Ui = x Uj to conclude that a Ui = a Uj . Hence, the self¬ 
connection functions in V 2 are all equal, and so in this case f{x) = J2 ee £ d( x p(e)> £ r (e))> f° r 
some smooth function 5. 

If Q is not the star graph, then there exists u G Vo with d(u) > 2. Consider then two 
edges (u, Vi), (u, V 2 ) G £. We have df/dx u invariant by interchanging v\ and V 2 , so 

dPv 1 / \ . dp v2 

which yields /3„ 1 = P V2 . Since £/ is connected, we apply this idea transitively through all 
edges of Q to conclude that 


x 


dp, 


dp v 


V 2 


■>Xu) (x V 2 ^xp) T (x Vll x u P 


dy 


dy 


f ( x ) = y y p(x v ,x u ) + y ^(^). 

ueVi ue ^( 1 ’) 
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Now let Vi,i ’2 G Vi with d(v i) = d{v 2 ). For any a,b, take 


x = (x Vl ,... ,x Vn ), x Vl =x V2 = a, u = 6, Wu G /(ui) U I[v 2 ). (2.6) 


By (11.41) applied to df /dx Vl and df /dx V2 at x, we have 

d ( v i)^( a ,b) + a' vi (a) = d(v 2 )^(a,b) + a' V2 (a), 
so a Vl = a V2 . Analogously, for v\,v 2 G V 2 with d(v 1 ) = d(y 2 ) we obtain 

d(vi)^(b, a ) +a' vi (a ) = d(v 2 )^(b, a) + a' V2 (a), 


and so a vi = a V2 . Therefore, self-coupling functions are invariant by degree on each subset 
of the partition of V, yielding (12.4p . 

Suppose now that there exist v\ G Vi and v 2 G Vo with d(y 1 ) = d{v 2 ) = d. Apply (jl.4[) 
to df /dx Vl and df / dx V2 at x given in (j2.6[) to obtain 


<9/3 d/3 

d — (a, 6) + a vii a ) = d—(6,a) + a v 2 \ 0 )i 


that is, 

d(^(a,b) -^(6, a)) = (a V2 - a Vl )\a). 

Since /3 does not depend on each variable independently, the above equality implies that (3 
is Z 2 -invariant and a vi = a V2 , resulting in (12.51) . 

Assume now that Q is not bipartite. We consider a spanning tree T of Q to rewrite (12.31) 
as _ _ 

/(•^) ^ ^ /9e(-^p(e)) *^r(e)) T ^ ^ /3e (^p(e)> ^r(e)) T 

ee£(T) e&£\£(T) veV 

Adding one edge e of ty in T will create a cycle, called a fundamental cycle. In fact, 
there is a one-to-one correspondence between edges e outside the spanning tree and the 
fundamental cycles C e of Q (see M). So we can rewrite the expression above considering 
edges corresponding to even cycles and odd cycles, that is, cycles with even or odd number 
of edges: 

f( x ) = /de(%p(e)i *^r(e)) T ^ ^ /3e(^p(e); ■^'r(e)) A ^ ^ Oi v (^hi); (2-7) 

ee£(H) ed£ veV 

where V. is the maximal subgraph of Q containing no odd cycles , T C % (f Q, and 


£ = {e G £\£(T) : Ce is an odd cycle}. 


T-L is bipartite, so it defines a partition V = Vi U V 2 so that 


f(x) 


^ ^ ^ ^ “1“ ^ ^ /^e(^p(e)? ^r(e)) “I - ^ ^ (^t?)• 

»eVi {u,v)^£(W) e^£ ve]3 


( 2 . 8 ) 
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Let e = (vi,V 2 ) G £. Notice that v { , v 2 G Vi or vi,v 2 G V 2 , otherwise Cg would be an even 
cycle. Consider (m^^i), (u 2 ,v 2 ) G £{T) neighbour edges of e. If v\,v 2 G Vi, the invariance 
of df/dx Vl by permutation of x Ul and x V2 and of df/dx V2 by permutation of x U2 and x Vl 
implies that 


df3 dfie 

dx ^ Xvi ’ Xui + dx ^ Xvi ’ Xv2 

dP ( \, d ^( \ 

dx V 2 ? ^ U 2 j %V2 / 


df3 dfe ( x 

9a; ' Et ’ 1 ’ Xv2 ^ 9a; ’ ‘ Eni ’ 

(x V2 , a; vi ) T \%U 2 > %V 2 ) • 


Hence, for any a, b, 


dfk 

dx 


(< a,b ) 


dp 

dx 


(a, b), 


dfe 

dy 


(b,a) 


dp 

dx 


(< a,b ). 


Therefore, /3g is Z 2 -invariant and /3g = p. If v\,v 2 G Vo, the result follows analogously. Finally, 
if ui,u 2 G V and d(v 1 ) = d(u 2 ), then at x as given in (12.61) we have df /dx vi (x) = df/dx V2 (x), 
and using the Z 2 -invariance of P we get 


dp dp 

diyP — ^.b) + a' Vl { a ) = d(v 2 )—(a,b)+a V2 (a), 


therefore a vi = a V2 yielding (12.51) . 


As a consequence of this theorem, a necessary condition for a network function / : P —> M 
to be admissible is 

^ ^ = 0, if (u,v) p £, (2.9) 


dx u dx v 


which captures the ‘dependency’ relation implied by the network structure (see m). 


Corollary 2.5. If Q is a regular graph, then the general form of an admissible function 
reduces to 

f( x ) = ^2P( x p(e), x r(e)), 

e^£ 

for some smooth Z 2 -invariant function <f> : K 2 ^ —» M. 

Proof: If Q is a regular graph, then self-connection functions are all identical. 

0 

From now on, whenever we label vertices with numbers, their variables shall be indexed 
with the corresponding numbers. However the subscript of a self-coupling function still 
denotes the degree of the corresponding vertex. 

We illustrate Theorem 12.41 with three graphs: 


Example 2.6. The admissible functions associated with each of the graphs presented in 
Figure [Til are given, respectively, as follows: 














1. fi(x) = rj(x 2 , Xi) + r](x 2 , x 3 ) + r](x 2 , £ 4 ), for some smooth 77 : R 2fc —» R. 

2. / 2 (x) = P(xi,x 2 ) + P(x 2 ,x 3 ) + /3(x 2 ,x 4 ) + /3(x 3 ,x 4 ) + a 3 (x 2 ) + a 2 (x 3 ) + a 2 (x 4 ), for 
smooth aj : R fe —> R, i = 2, 3 and (5 : R 2fc —» R Z 2 -invariant . 

3. / 3 (x) = 0(xi, x 2 ) + 0(xi, x 5 ) + 0(xi, X 6 ) + 0(x 2 , x 3 ) + 4>{x 2 , x 4 ) + (fi{x 3 , X 4 ) + 0(x 3 , x 5 ) + 
0 (x 4 ,X 6 ) + 4>(x5,Xg), with <p : R 2fc —>• R Z 2 -invariant . 




Qi 


$2 



£3 


Figure 2.1: Graph t?, corresponding to the admissible function /j given in Example 12.61 z = 1,2, 3. 

We just notice that we first use Theorem 12.41 to write f\ as 

fi(x) = f3(x 2l xi) + /3 (x 2 ,x 3 ) + fi(x 2 ,x 4 ) + ai(x 2 ) + «i(xi) + ou(x 3 ) + ol\{x 4 ), 
and then define 77 ( 0 , 6 ) = @(a,b) + aa( 6 ). 

We give below another example of a bipartite graph, with a larger number of cells: 

Example 2.7. We consider a bipartite graph for which d(u) 7 ^ d(v) if u and v have distinct 
colours, see Figure [2721 So an admissible function for this graph is of the form (12.4ft and it 
is given by 

/(*) = 5Z P{x 2 ,Xi)+ ^2 P{x4,Xj)+ ^2 P(x 6 ,X k ) +a 1 (x 1 ) +ai(xi 0 ) + 

i=l,3,5,7,8,9 j= 3,5,7,8,9 fe=3,8,9,10 

+ a 2 (x 5 ) + a 2 (x 7 ) + a 3 (x 3 ) + a 3 (x 8 ) + a 3 (x g ) + 74 (x 6 ) + 7 5 (x 4 ) + 7 e(^ 2 )- 
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3 Critical points on regular graphs 

In this section we direct attention to regular graphs, so admissible functions are of the form 

f( x ) = x r(e)), (3-1) 

ee£ 

with 0 : M? k —» R Z 2 -invariant. Our aim is to discuss when the nature of certain types 
of critical points of / can be deduced from the associated critical points of the coupling 
function 0. Totally synchronous critical points and 2-colour patterns of critical points are 
considered, and these are given in Subsection 13. II and in Subsection 13.21 respectively. Minima 
of one function are detected from minima of the other depending on the graph architecture. 
It is clear that our analysis is distinct for these two types of critical points in the class of 
admissible functions. Critical configurations with three or more colours must be discussed 
under a different approach, using results in [lj. 


3.1 Totally synchronous patterns 

For an admissible function of the form (13.11) . a point x° = ( x 0 ,..., x 0 ) is a critical point of / 
if, and only if, (x 0 ,a;o) is a critical point of its coupling 0, even when 0 has no permutation 
invariance. In this subsection we assume that the graph is regular and that the configuration 
space of each cell is one-dimensional {k = 1). In this case, the coupling is Z 2 -invariant, and 
so the nature of such critical points given by the Hessian criterion is then a 2-parameter 
problem, the parameters being given by the second partial derivatives of 0 at (xq^o) G R 2 . 


We shall use below subscripts 1 and 2 for the partial derivatives of 0 with respect to 
the first and second variable respectively. In what follows we denote a = 0n(a;o,a;o) and 
/3 = 0i 2 (x o ,xo). The Hessian of / at x° is the matrix H = (a + /3)dl — (3£, where I is the 
identity matrix and C = D + A is the graph Laplacian, where D is the degree matrix and A 
is the adjacency matrix of Q ([2D]). Notice that the eigenvalues of H are 

Hi = (a + /3)d - (3\i, (3.2) 


where 0 = Ai < A 2 < ... < = A are the eigenvalues of C. 


Proposition 3.1. Let Q be a d-regular graph and consider an admissible function on Q of 
the form H3.1\ ) . Then generically a totally synchronous critical point x° = (x 0 ,... ,x 0 ) is a 
minimum of f if, and only if, 

0 n(a? 0 ,ar 0 ) + (fi 2 (x 0 ,x 0 ) > 0 , <i>u.{x 0 ,x Q ) + (1 - 3 ) 012 (^ 0 ^ 0 ) > 0 , 

d 

where A is the greatest eigenvalue of the Laplacian of Q. 


Proof: The Hessian H is generically nondegenerate, in which case x° is a minimum if, and 
only if, 


a > 


—d + A i 
d 


P, 
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Figure 3.1: The wedge region where the converse of Corollary 13.21 fails. 

for all i — 1,..., n. It is well-known that A < 2d, where equality holds if, and only if, Q is 
bipartite. Hence, for any i = 1,... ,n, we have (—d + Xfi/d e [—1,1], and so a minimum 
occurs if, and only if, 

a + (3 > 0 and a + (—— -)/3 > 0. (3.3) 

d 

□ 

Inequalities (|3.3[) are precisely the Hessian criterion for a minimum of 0 in (13.11) if A = 2d. 
So, for the generic case when this criterion applies, we have: 

Corollary 3.2. An admissible function f defined on a regular graph, given by h3.1\) . has a 
minimum at x° = (x 0 , ... ,£ 0 ) if the coupling function 0 has a minimum at {;xq,Xq). The 
converse holds if, and only if, Q is bipartite. 

Looking at the plane of the parameters (pn(x 0 ,x 0 ), (pi 2 {x 0 ,x 0 ), the converse in the corol¬ 
lary above fails precisely in the wedge 

011 Ob), xq) - 0120CO, x 0 ) < 0 , 4>n(x 0 , x 0 ) + (1 - 4)012(^0, x 0 ) > 0 , (3.4) 

d 

see Figure [3711 Now, for a d-regular graph, we have the lower bound for the greatest eigen¬ 
value A, 

A > d+ 1, 

where equality holds if, and only if, the graph is complete [23]. This means that we can 
measure the set in parameter space where minima of / are not in one-to-one correspondence 
with minima of the coupling function 0. Its ‘size’ is determined by the topology of the graph: 

Corollary 3.3. For regular graphs, on the plane of the two parameters (pii(xo,x 0 ) and 
0i 2 (xo,xo) the wedge determined by vanishes if, and only if, the graph is bipartite. 

It is as large as possible if, and only if, the graph is complete. 

More generally, for any regular graph we can use (13.21) to track down the nature of a fully 
synchronous critical point x° of / for different values of a and (3. Equivalently, we detect 
the degree of instability of x° thought of as an equilibrium of the vector held —V/. In fact, 
in the open region (a + f3)(a — fi) > 0, x° is a minimum of / if a > 0 and it is maximum if 
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a < 0. If n is the number of cells in Q , we have that in the region (a + /3)(a — ft) < 0 are 
the lines on which fij = 0 , for 1 < j < n, whose equations are 


with algebraic multiplicity of fij given by the algebraic multiplicity of X 3 for L. On this line 
and for j3 > 0, fik > 0 if 1 < k < j and fik < 0 if j < k < n; for j3 < 0, ptk < 0 if 1 < k < j 
and Hk > 0 if j < k < n. Also, the line a + /3 = 0 corresponds to the case for which H is 
a (3 -weighted Laplacian, that is, a weighted Laplacian (see [20]) with all weights equal to f3. 
So Hi = 0 and a minimum of / can only occur for negative values of /3. Finally, on the line 
a — /3 = 0, for positive (negative) values of a we have /q- > 0(< 0), for 1 < j < n unless the 
graph is bipartite, in whose case all eigenvalues are still positive (negative) except fx n , which 
vanishes. 

With the two corollaries above in mind, it becomes natural to look for a class of nonregular 
graphs for which the wedge vanishes, that is, for which the converse of Corolllary 13.21 does 
hold. As we see in the next proposition, this is the case for the complete bipartite graph 
K m n , i.e., the bipartite graph where every vertex of one set of vertices Vi is attached to 
every vertex of the other set V 2 , with |Vi| = m, | V 2 1 = n. If m — n, K myn is regular and it 
fits in Corollary 13.21 so in the next result we are interested in the case m 7 ^ n. 

Proposition 3.4. For distinct m,n > 2, an admissible function f given in h3.1\) defined on 
the complete bipartite graph K m:Tl has a minimum at x° = (xo,... ,xo) if, and only if, the 
coupling f has a minimum at {xq,xq). 

Proof: The Hessian matrix of / at x° is given by 



where a = 4>i\{xq,xq), /3 = <f 12 (^ 05 ^ 0 ) and B is the n x m-matrix in the adjacency matrix 



(3.5) 


of K mn , with all entries equal to 1. We shall now describe the elements in the spectrum 
Spec (H(x 0 )). For A 7 ^ na, we have that det(//(x°) — AI m+n ) = 0 if, and only if, 



The matrix in the expression above can be rewritten as 

I n A 2 — (m + n)al n X + mna 2 \ n — /3 2 BB 1 = (M + — AI n )(M_ — AI n ), 


where 
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which are matrices with real entries since BB t is a positive semi definite matrix. We now 
investigate the eigenvalues of these two matrices distinct from na. We notice that the 
eigenvalues of BB f are given by p, = ( 2 , for each eigenvalue f of the adjacency matrix A. 
But it is well-known that ( = 0, — y/mn, \Jmn. Hence, the eigenvalues of M + and M_ are 

1 / ,\l/2 

- [( m + n)a ± ^(m — n) 2 a 2 + 4 p/3 2 J ], 


for each of the two values of p. Therefore, the eigenvalues distinct from na are ma and 


1 

2 


[(m + n)a ± 



For A ma we proceed analogously as above for A 7 ^ na to conclude that na is also an 
eigenvalue of H(x°). And now it is straightforward to check that for any m and n these 
eigenvalues are all positive if, and only if, (a + (3), (a — (3) > 0. g 


3.2 2-colour patterns 

In this subsection we investigate critical points x of admissible functions (13. ip on regular cell 
graphs Q which are not totally synchronous and which come from critical points ( Xo,yo ) of 
the coupling function with x 0 7 ^ yo■ So, for V = { 1 , 2 , ..., n}, we are considering 

x=(x!Xi = Xq or ?/o, i = l,...,n. (3.6) 


We start by showing that no such critical points are expected unless the graph is bipartite, 
even if it is non-regular: 

Proposition 3.5. Let Q = (V,£) be a cell graph. For a generic coupling function </> a point 
x G M fcn as in H3.6\) can be a critical point of an admissible function of the form 

f( X ) = (3.7) 

e^S 


<f : W 2k —> R, only if Q is bipartite. In this case, for V = Vi U V 2 , x is given by 


Xi = 


Xq, i G Vi 

yo, i G V 2 . 


(3.8) 


Proof: We may assume that the coupling function </> is Z 2 -invariant; otherwise, the graph 
is already bipartite, by Theorem 12.41 Suppose that (x 0 ,y 0 ) is a critical point of <f with the 
generic condition that neither are ( Xq,Xq ) nor (yo,yo). Fix j, 1 < j < n. If x 3 = x 0 , let p be 
the number of vertices £ G I(j) such that X£ = Xq and q be the number of vertices £ G I(j) 
such that xi = yo. Then for any i, 1 < i < k, 

,X ^ = P ^( Xo ’ x °) + q ^i( x o, y o) =P ^-(xo,xo). 
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By the permutation invariance of 0, <90/&c* (x 0 , x 0 ) 0 for some i, 1 < i < k. Hence, p — 0, 

and therefore X£ = r/ 0 for all l G J(j). Analogously, if Xj = y 0 , then X£ = x 0 for all £ G J(j). 
This dehnes a partition of V so that Q is bipartite and the components of x must be given 
by (USD. D 

As a consequence of the proposition above, the assignment of one distinct colour for each 
x 0 and y 0 provides a standard 2-colouring on the graph, that is, no two vertices sharing 
the same edge have the same colour. Let us also notice that in the case when the coupling 
function is invariant by permutation, then 2-colour critical points of / come in pairs. 

If Q is a regular bipartite graph, then it has an even number of vertices, n = 2m, m > 1, 
with |Vi| = |V 2 1 = m. I 11 fact, just impose the regularity condition on the well-known 
equality for bipartite graphs, 

Y de s( w ) = de s( u )- 

ugV 1 ugIA 

So here we shall investigate 2-colour patterns constrained to the set of (d, m)-graphs, as we 
define next: 

Definition 3.6. A connected graph Q = ( V,£ ) is a ( d,m)-graph if it is a d-regular bipartite 
graph with |Vi| = |V 21 = m- 

We refer to HIE] for the number of existing (d, m)-graphs when d = 3,4, 5 and 3 < 
m < 16. Figure [3T2l illustrates the cubic graph, which is a (3,4)-graph. 



Figure 3.2: The cubic graph is a (3,4)-graph. 


I 11 the next two propositions we assume the variable on each cell is one-dimensional. 


Proposition 3.7. The spectrum of a 2-colour critical point of an admissible function depends 
only on the coupling function and the integers d and m. 


Proof: Given any (d, m)-graph Q = (V, £), V = Vi U V 2 , use {1,2,..., m} to label elements 
in Vi and {m + 1, m + 2,..., 2m} to label elements in V 2 . For any admissible function (13.ip 
on Q , its Hessian at the 2-colour critical point x given by (13.81) is of the form 


where a = (fii(x 0 , y 0 ), P 
which is of the form 


H{x) = 


A = 


dal m 

0 

0 

dqlm 

— 011 (do, %o) 

~ f ° 

C l \ 


0 ) 


+ /3A, 


(3.9) 


0i2(•£(); do) j 7 = 0ii (yo, xq) and A is the adjacency matrix of Q, 


(3.10) 
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for an order-m matrix C. 

Any other (d,m)~ graph Q is given by permuting edges of Q, so there exists an order-m 
permutation matrix P such that its adjacency matrix is given by 


A = 

But the Hessian matrix for Q at x is 

H(x) = 


( 0 

(.pcy 

y PC 

0 

dotl-rn 

0 

0 

d'ylm 


+ ft A, 


which is similar to f!3.9f) : just take 


as a similarity matrix. 


P 

0 

0 

Ira 


0 


The next result is the analogous of Proposition 13.41 for (d, m)-graphs at 2-colour pattern 
configurations. 


Proposition 3.8. If Q is a ( d,m)-graph, then a 2-colour pattern x in Q as in H3.8\) is a 
minimum configuration of an admissible function given in \ 3. 1\) defined on Q if, and only if, 
(xo,yo) is a minimum of its coupling function. 


Proof: The Hessian matrix of / at x is given by fl3.9p - fl3.10p . and its characteristic polynomial 
Ph{x) is given by the product of two characteristic polynomials 


Ph(x) — Pn + Pn_, 

where 

1 / 1/2 
N± = - [d(a + 7 )I m ± [d 2 {a - 7 ) 2 I m + ApCC') ], 

which are matrices with real entries since CC* is a positive semidehnite matrix. Hence, 

Spec (H(x)) = Spec(Al + ) U Spec(A^_). 

Now the least element of Spec (H(x)) is the minimum of Spec(Al_), which is given by 

^ (d(a + 7 ) - ( d 2 (a - 7 ) 2 + Aufi 2 ^j 1 ), 

where v is the greatest eigenvalue of CC l . But v = A 2 , where A is the largest eigenvalue 
of the adjacency matrix of Q. Since Q is d-regular and bipartite, it follows that A — d. 
Therefore, 

6v_ = ^ ((a + 7 ) - ((« - 7) 2 + 4/3 2 ) 7 ), 

which is positive if, and only if, a > 0 and cry — /3 2 > 0. g 

We end this section with a remark regarding graphs in general, not necessarily regular. 
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Remark 3.9. As a consequence of Proposition 13.51 unless there are additional symmetries 
other than possible permutation invariance of the coupling function, occurrence of both 2- 
colour patterns and totally synchronous critical points (zq, , z 0 ) is expected for admissible 
functions (13.7ft only if zq ^ Xo,yo. However, under extra symmetry constraints, these two 
types of configurations can co-exist for zo = xo or y 0 . In fact, these can be critical config¬ 
urations, among others, for example when the admissible function has an invariance under 
the circle group S 1 . In this case these two can even both be minima of /. We shall return 
to this point in the next section. A study of a particular SMnvariant function on a ring is 
carried out in Subsection 14.21 

4 S^invariant admissible functions 

In this section we analyse possible critical configurations of a class of admissible functions 
for which the interest lies in the behaviour of phase differences between planar unit vector 
states of coupled cells. For each cell i in the graph we can identify the vector state with 
its angle 6 i with respect to a fixed direction, which we take to be the vertical direction. In 
this setting, for a network represented by a graph Q of n cells, the configuration manifold is 
the n-torus T n = S 1 x • • • x S 1 . In this case the admissible functions on JV = ( Q,T n ) are 
additionally invariant under the translational diagonal action of S 1 , 

h{ 6 1 + • • • 0 n + ip) = h{9 1, ...,9 n ), £ S 1 . 

Imposing this condition on the general form given by Theorem 12.41 yields 

®d(v) = 0, Vv G V, (4.1) 

so an S 1 -invariant admissible function has only trivial self-connections. In addition, the 
coupling /3 in (12.41) or (12.51) is of the form 

P(x,y) = S(x-y), (4.2) 

for some smooth function <5 : S 1 —> R. If 5 admits an extension to M. then this must be 
an even function for the case (12.51) . We remark that (14.ip is a consequence only of the S 1 - 
invariance, so it also holds for SMnvariant functions defined on non-regular cell graphs (see 
Corollary 12.5[) . 

Functions of this type appear in many applications, including the Kuramoto model (see 
|H [IB]) and the Antiferromagnetic XY model (see [121 EH H])• We briefly describe the 
second, which models alignments of spins in magnetic materials. The spins are distributed 
in a planar lattice, the ‘site’ of each being each vertex of the lattice. The source of the general 
behaviour of magnetic materials is the spontaneous, parallel alignment of neighbouring spins. 
This type of magnetism is called ferromagnetism. An equally common yet very different form 
of magnetism is antiferromagnetism. Like ferromagnetism, this arises due to the spontaneous 
alignment of neighbouring spins, but the alignment of spins in antiferromagnetic materials 
is anti-parallel (in an up-down configuration). The total energy of this spin system, the 
Hamiltonian, assumes a particular form of an admissible function defined on the lattice, 
with coupling functions typically being given by the cosine function and the sum related to 
couplings is over nearest neighbour spin sites. 
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4.1 The Hessian 

Let us first consider the case of nearest-neighbour coupling of n identical cells in a ring, so 
the cell graph Q is a regular n-sided polygon. The admissible function is then given by 

n 

h(9) = (4-3) 

i =1 

where 9 0 = 6 n and 5 is even. 

Critical points of h are points 9 — (9 1 ,..., 9 n ) such that 

8\9 i+ 1 - 6i) = S'(9 i - 9^), l<i<n-l. (4.4) 

The Hessian matrix H of h has its entries given as follows: for 1 < % < n, 

5 = 6"(9 l+1 -9 l ) + 6"(9 i -9 i _ 1 ), 

au.i 

d 2 h = f -5%9 i - 9 ,), j = i + 1 , i - 1 

09,09-j \ 0, otherwise. 

From this particular case of a ring it is straightforward to see how the expression of the 
Hessian generalizes for any cell graph Q of n vertices: If Q is bipartite (and 5 may not be 
even), from the partition V = Vi U V 2 of the set of vertices the admissible function is of the 
form 

h(9 l ,...,9 n ) = J2 (4.5) 

ieV fe , jei(i) 

for arbitrary k — 1 or 2. Set l = 1 or 2, l ^ k. Then the critical points 9 of h have coordinates 
that solve the equations 


X s '(»< 

~0j) = 

0, 

for % 

£ Vfc 

ie/(i) 





£ Wi 

-9i) = 

0, 

for i 

6 V,. 

jei{i) 





as follows: 

for 1 < i 

< 

n, 


h _ J £* 

E/(0 

— 

9j ), if i e 

f “ l £* 

•-1 

— 

9i), if i G 

r 

-Oj), 

if 

i e V fc , 

j e 

{ -<m 

- 0i), 

if 

ieVi, 

j e 

1 0 , 


otherwise. 



09? 


0 2 h 

09,09 j 

If Q is not bipartite, then the function is of the form 

h(9 u ...,9 n ) = Sidi-Oj), 


(4.6) 


(4.7) 


where 6 is even and the equations for critical points as well as the Hessian entries above 
simplify in the obvious way. 
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Example 4.1. Consider the graph in Figure 14.11 An admissible function defined on this 
graph is of the form (14.51) and the Hessian matrix is given by 


/ 

«21 

—«21 

0 

0 

0 

0 

0 

\ 


—«21 

021 + 023 

—023 

0 

0 

0 

0 



0 

—023 

023 + O43 + CJ53 

—O43 

-053 

0 

0 



0 

0 

—043 

043 + “46 

0 

—046 

0 



0 

0 

—053 

0 

053 + a 57 

0 

—057 



0 

0 

0 

— O46 

0 

046 

0 


V 

0 

0 

0 

0 

—057 

0 

057 

/ 


where = 5"{0i — Of). 




Figure 4.1: A bipartite graph and its 2-colouring. 

Remark 4.2. For any cell graph it follows from the description above that at any critical 
point the Hessian of an admissible SMnvariant function h is a weighted Laplacian. When 5 
is an even function, the weight at each (i, j)-edge is 5"(6i — 9f. If 5 is not even, then Q is 
bipartite and h is given by (I4.5jl : in this case the weight at the (i, j)-edge is still 5" {9, — 9f), 
but now i G 14. 

It is well known that without symmetry constraints, a smooth function on a manifold 
generically is a Morse function. In the present case, the Hessian, being a weighted Laplacian, 
has a zero eigenvalue (with eigenvector with all elements equal to 1). We then have: 

Proposition 4.3. For any cell graph Q, all critical points of S 1 -invariant admissible func¬ 
tions on Q are degenerate, so no Morse singularities are to be expected. 

Let us now remark that simultaneous existence of totally synchronous and 2-colour critical 
points can be expected due to extra symmetries; see Remark 13.91 In view of Proposition 13.51 
in the search of 2-colour critical points bipartiteness of Q is a generic assumption. In this 
case, V = Vi U V2 and 2-colour critical points are given by 1 9 — {9 1,..., 9 n ), 

n _ f a 0> R V 1 

‘ “ I An V 2 , 

a 0 f (5 0 . Then, we have synchrony and a 2-colour configuration simultaneously as long as 
the coupling function has two distinct critical points, one of them being the origin. I 11 fact, 
9 0 = (9 0 ,..., 9f) is critical for h if, and only if, 5'(0) = 0 and, from (14.61) . 9 is critical for 
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h if, and only if, 8'(a o — /3q) = 0. In addition, the Hessians at 9° and at 6 are the weighed 
Laplacians with weights 5"(0) and 8"(ao — /do); that is, if C is the graph Laplacian of £/, then 

H(6°) = 8"{ 0)£, H(6) = 8"(a 0 -(3 0 )£, 

and hence the nature of one pair is determined by the nature of the other. The aim of the 
next subsection is to consider the case when Q is a ring and the coupling function is under 
some special constraints. We shall encounter further critical points co-existing with these 
two, and their nature shall also be investigated. 

4.2 Stable synchronous patterns on a ring of many cells 

We assume that the graph is a ring of n cells, so the admissible function is given by ( 14 . 3 [) : 

n 

h(0) = J2 s &- e i- 1 ), 

i=1 

with #o = 0 n . By rescaling, we consider here S 1 = [0,1] and the domain of 8 in (14.2p to be 
the unit interval, 8 : [0,1] —» M. Assume that 8 can be extended to an even smooth function 
to M. On the configuration space R n , the gradient of h is an admissible vector field on the 
n-ring. In applications, h represents the total energy associated to the gradient system 

6 = -V/i(0), (4.8) 

whose stable equilibria correspond to minimum points of h. Based upon that, the global 
minimum of h shall be called here the ground state. From Proposition 14.31 the critical points 
of h are not Morse singularities, for its Hessian at any point is a weighted Laplacian, and so 
the kernel is at least one dimensional. With that in mind, we shall say that an equilibrium 
of (14. 8p is asymptotically stable , or simply stable , if all nonzero eigenvalues of the Hessian of 
h at this point are positive. 

As mentioned previously, in many applications the coupling is the cosine function; see 
Here we consider a slightly more general class of coupling, but chosen 
suitably to provide a clear understanding of how it combines with the architecture of the 
graph to reveal the ground state for the system other than the obvious alternating up-down 
configuration, which corresponds to antiferromagnetism interactions mentioned earlier. More 
precisely, we assume that the coupling function 8 satisfies the following conditions: 

(Cl) 8 is even and 1-periodic; 

(C2) The first derivative 8' vanishes only at 0 and 1/2; 

(C3) The second derivative 8" is monotone in the interval (0,1/2). 

We now introduce the variables 


Xi = 9 i+ 1 -9i, 1 < i < n, 
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where 6 n+ 1 = 0\. From (Cl) we can restrict our analysis to x* G [0,1/2], for i — 1,..., n. 
Critical points of h are then given by 8'(x i+ 1 ) = 8'{xi ), 1 < i < n — 1, see (14.411 . where 

n 

Xi G Z. 

Z=1 

For m denoting the sum above it follows that 0 < m < n/2. In addition, due to the S 1 - 
invariance of h it suffices to describe its critical points on T n . In fact, for the representation 
of the dihedral group D ra of permutations that preserve the n-ring, this is a D n x SMnvariant 
function, for the D n x S 1 -action on T n given by 

(7T, p) • ($ 1 , • • • , d n ) i(i) “I - ip , • • • , 0 n — !(n) 4“ P) • 

In Table Q] we present the list of all possible equilibria up to their isotropies in D n x S 1 . We 
denote by Z n (r, ip) the isotropy subgroup generated by (r, ip), where r is the permutation 

r = f 1 2 "• " 

V 2 3 ■" 1 

and p = 7n/n, m = 0,..., \_n/2\. If we label the cells clockwise, a Z n (r, (^)-isotropy equilib¬ 
rium is a configuration such that each cell has a phase shift by p with respect to its prior 
neighbour (Figure l4~2l on the left), whereas at a trivial isotropy equilibrium each cell has a 
phase shift by £ or r] relatively to its prior neighbour, where £ ^ 7/ G [0,1/2]. As we shall 
see in Proposition 14.41 below, a trivial isotropy equilibrium is stable only if its configuration 
is of the form given in Figure l4~2l on the right for £ > r]. 


6 


n 






o 2 


<p 




£ 


On-1 


O3 0 n —i 


03 


Figure 4.2: On the left, a Z n (r, </?)-isotropy equilibrium, for which 6i + \ = 6i + p, i = 1,... ,n with 
9 n+ 1 = 9 1 . On the right, a 1-isotropy equilibrium, for which 9i+\ = 9i + £, i = 1,..., n — 1 and 

9i = 9 n + V, 9 ^ £. 

It is a direct consequence of the conditions (C1)-(C3) that there are no equilibria rather 
than those in Table Q] Also, we notice that trivial isotropy equilibrium points only occur for 
special couplings. In fact, the values £ and 7] in the last row of this table are such that 

0 < V < £ < 1/2, 5\0 = 8'(ri), (4.9) 
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whose existence is guaranteed by the conditions imposed on 5. However, these will correspond 
to equilibria of the system if and only if p£+qr) is an integer m, for p, q and m in the conditions 
presented therein. This is the case for example for the values £ = 1/2, r) — 0, giving rise to 
equilibrium points 6 = ( 61 ,, 6 n ) such that there is en even number of pairs of neighbours 
whose states differ by 1/2. 


m € Z, 

0 < m < n/2 

Isotropy 

Critical point 
representative 

Value of 

energy 

Laplacian 

weights 

0 

D n 

o 

o' 

o' 

n 6(0) 

v'(o) 

^ 0 , n/2 

Z n (r,m/n ) 

(0, m/n, 2m/n, ..., m) 

n 6(m/n ) 

6" (m/n) 

n/2 (if n even) 

Z„(r,l/2) 

(0,1/2,0,..., 1/2) 

n 6(1/2) 

6"( 1/2) 

+ qq, 

for p, q G N, p + q = n 

1 

(#1) ■ • • ,9 n ) 

9i+ 1 - 0i = £ or 7] 

p6(£)+q6(rj) 

s"(v) 


Table 1: Critical points of admissible functions on a ring of n cells. 


By direct investigation, if ^(O) > 0 then the fully synchronous pattern Xi — 0, 1 < i < n, 
is the ground state. This corresponds to the ferromagnetic state and it falls into the case 
a + /3 = 0,/3<0 discussed in Subsection 13.11 

On the other hand, if ^'(O) < 0 then fully synchrony corresponds to a totally unstable 
equilibrium. The aim here is to extract from Table [Q all possible stable equilibria in this 
case. In addition, we detect the ground state if we allow the ring to have a large number of 
cells. In this context, we shall then consider: 

(C4) S"(0) < 0 and S"(l/2) > 0. 

The results are given in Propositions 14.41 and 14.61 and rely on the fact that the Hessian 
of admissible functions defined on rings is, at any point, a graph weighted Laplacian. We 
shall see how stability is attained depending on the topoloy of the graph combined with its 
weights. The weights are real numbers associated with each edge, so as expected they are 
determined by the constrains imposed on the coupling function 5. 

We first recall an important fact of [5j about eigenvalues of weighted Laplacians. For 
any weighted graph Q with n vertices, let Q + (resp. Q_) denote the subgraph with the same 
vertex set as Q together with the edges of positive (resp. negative) weights. We also define 
the three indices no,n_ and n + as the numbers of zero, negative and positive eigenvalues 
of the weighted Laplacian. In [5] the authors give the best possible bounds on these three 
indices involving only topological information, i.e., connectivity of the graph and the sign 
information on the edge weights. For c(Q) denoting the number of connected component of 
any graph Q , the estimates are given as follows: 

c(G+) — 1 < n_ <n — c(G~) 

c(G -) — 1 < n + < n — c{g+) (4-10) 

1< n 0 < n + 2 - c(G^) - c(G + ). 
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Proposition 4.4. Let Q be a ring of n cells coupled under conditions (Cl)-(Cf). Then the 
stability of critical points that appear in Table{J\ is as follows: 

(1) Forn even, an equilibrium with Z n (r, 1/2) -isotropy is stable; 

(2) Form / 0, n/2, an equilibrium with Z n (r, m/n)-isotropy is stable if and only if 8"(m/n) > 

0 ; 

(3) For n odd, a trivial isotropy equilibrium 6 = (0 1; ..., 0 n ), 9., e {0,1/2}, is stable only if 
at most one pair of neighbours assume the same state value. 

(4) Possible trivial isotropy equilibria with f,r] in the open interval (0,1/2) are stable only if 
p = n — 1, q — 1. 

Proof: The weighted Laplacian at the Z n (r, l/2)-isotropy equilibrium (for n even) is given by 
5"{1/2)C, where £ is the graph Laplacian (without weights). The n— 1 nonzero eigenvalues 
of L are positive, so stability of this point follows from (C4). The stability of Z n (r, m/ri)- 
isotropy points follows analogously. Consider now the trivial isotropy points 9 given in 
(3). The weights in this case are given by 5 7/ (0) and 5"{ 1/2). Now, from (14.101) a necessary 
condition for any point to be stable is c(Q + ) = 1, i.e., Q + must be connected. Under condition 
(C4) this is the case on a ring if and only if the states of cells alternate between 0 and 1/2 
except for possibly one pair of neighbours. Finally we consider 9 given in (4). We have that 
5"(£) > 0 and 8"{rj) < 0 by (C3), and the result follows analogously as for case (3). g 

Remark 4.5. Case (3) of Proposition 14.41 corresponds to the frustration phenomenon in 
magnetic systems. Geometrical frustration is a feature in magnetism related to the arrange¬ 
ment of spins. In a ring with an even number of spins, the energy is minimized when each 
spin is aligned opposite to neighbours (see Proposition 14.61 below). If the number of cells is 
three for example, once the first two spins align anti-parallel, the third one is ‘frustrated’ 
because of its two possible orientations, up and down, which give the same energy. The 
third spin cannot simultaneously minimize its interactions with both of the other two. From 
Proposition 14.41 (3). this is the case if there is exactly one pair of neighbours with same state 
value, as illustrated in Figure EOl 


t 



4 


4 — t 


Figure 4.3: The frustration phenomenon in a ring of five cells. 

Proposition 4.6. Let Q be a ring of n cells coupled under conditions (Cl)-(Cf). If n is 
even then the ground state (9 1 ,..., 9 n ) of (| f.8\ ) is given by 9 l+ i — 9i — 1/2, i — 1,..., n. If 
n is odd and sufficiently large, then the ground state is given by 9 l+i — 9 t = 1/2 — 1/2 n, 
i — 1 ,..., n. 

Proof: For n even, the conclusion follows by direct investigation on the values of the energy 
function h given in Table |T| For n odd we have to compare among Z n (r, m/n)-isotropy and 
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trivial isotropy stable equilibria given by Proposition 14.41 Since 5 is a decreasing function on 
[0,1/2], the lowest energy with Z n (r, m/n) isotropry is attained for the largest m, namely 
m = (n — l)/2. This corresponds to the equilibrium 9 = (6>i,..., 9 n ), 9i = 1/2 — 1/277, 
i = 1 ,..., 77 . Hence, it remains to compare the energy values at this point and at trivial 
isotropy points 6 which are given by (3) and (4) in Proposition 14.41 

We have that (n — 1)£ + 77 = m, 0 < 777 < n/2 and 7 ? = £ — e, 0<e< 1/2. For n odd, 
777 < (n — l)/2. If 777 < (n — l)/2, then 

72 — 1 

(77 - l)f + 77 < —-1, 

so 

77 -I e - 1 1 1 

£ < -+-<- . 

2 7i 77 2 2n 

Hence, £(£) > 5(1/2 — 1 /( 277 )), and since 5(r]) > 0, 

< („-i)i( ? )+%). 


If 777 = (77 — l)/2, then 


77 —I e 77 — 1 1 — 77 

£ = + - V=~^ + -e- 


2n 77 


2/7 77 


If 77 is sufficiently large, then 


v 2 2n’ ’ v 2 2n n’ v 2 2n J 


We now just compare 




'2 2n / v 2 277 ; 


with 

f(9) = (77 - 1)5(0 + 5(7/) 

to conclude that, in this limit, h(6 ) < h{6). g 

Proposition 14.61 asserts that for the bipartite case the ground states are achieved in the 
up-down configuration, see Figure ICT ab This is the corresponding result for rings to the 
analysis in |13] of the classical AFXY model, for which bipartite lattices of spins attain 
the ground state when the two sublattices are aligned in opposite directions. For the non 
bipartite case, (when n is odd), the ground state value depends on n and tends to the 
up-down arrangement in the limit 77 —> 00 , as depicted in Figure l4~4T bh 
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0 1 — 03 — ■ ■ ■ 9-n-l 


02 ~ $4 — • • • — 02 n 


(a) 



Figure 4.4: Ground state in a ring of n cells under conditions (C1)-(C4) when (a) n is even and 
(b) n is odd and large. 
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